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ABSTRACT 

Aims. Heating associated with non-equilibrium nuclear reactions in accreting neutron-star crusts is reconsidered, taking 
into account suppression of neutrino losses demonstrated recently by Gupta et al. Two initial compositions of the nuclear 
burning ashes, Ai — 56 and Ai = 106, are considered. Dependence of the integrated crustal heating on uncertainties 
plaguing pycnonuclear reaction models is studied. 

Methods. One-component plasma approximation is used, with compressible liquid-drop model of Mackie and Baym to 
describe nuclei. Evolution of a crust shell is followed from 10^ g cm~^ to lO'^^'® g cm~^. 

Results. The integrated heating in the outer crust agrees nicely with results of self-considtent multicomponent plasma 
simulations of Gupta et al.; their results fall between our curves obtained for Ai = 56 and Ai = 106. Total crustal heat per 
one accreted nucleon ranges between 1.5 MeV/nucleon to 1.9 MeV/nucleon for Ai = 106 and Ai = 56, respectively. The 
value of Qtot depends weakly on the presence of pycnonuclear reactions at lO^'^ — lO'^^ gcm~^. Remarkable insensitivity 
of Qtot on the details of the distribution of nuclear processes in accreted crust is explained. 

Key words, dense matter - equation of state - stars: neutron - stars :general - X-rays: bursts - X-ray: binaries - nuclear 
reactions 



of many nuclides, and included the contribution from the 
nuclear excited states. They found that electron captures 
in the outer crust lead mostly to the excited states of the 
daughter nuclei, which then deexcite heating the matter. 
Consequently, they found that the neutrino losses in the 
outer crust were negligible, w hich strogly increased the 
outer crust heating, compared to lHaensel fc Zdunik I (|2003l ) 
(by a factor of four). However, outer crust contributes only 
a small fraction of the Qtot, and neutrino losses in the in- 
ner crust, where the bulk of Qtot is produces ar e smal l, 
so that the original value of iHaensel fc Zduniki ()1990[ ) , 
Qtot ~ 1-4 MeV/nucleon, remains quite a reasonable es- 
timate (see Sect. S] of the present paper). 

The phenomenon of deep crustal heating appears to be 
relevant for the understanding of the thermal radiation ob- 
served in the soft X-ray transients (SXTs) in quiescence, 
when the accretion from a disk formed of plasma flow- 
ing from the low-mass companion star is switched off or 
strongly suppressed. Typically, the quiescent emission is 
much higher than the expected one for an old neutron star. 
It has been suggested that this is due to the fact that the 
interiors of neutron stars in SXTs are heated-up, during 
relatively short periods of accretion and bursting, by the 
non-equilibrium processes associated with nuclear reactions 
taking place in the deep layers of the crust (Brown et al. 
1998). The deep crustal heating, combined with appropri- 
ate models of neutron-star atmosphere and interior, is used 
to explain observations of SXTs in quiescence. The lumi- 
nosity in quiescence depends on the structure of neutron- 
star core, and particularly on the rate of neutrino cool- 
ing. This opened a new possibility of exploring the internal 
structure and equation of state of neutron stars via con- 
frontation of theoretical models with observations of qui- 



1. Introduction 

Neutron star crust that is not in full thermodynamic 
equilibrium constitutes a reservoir of energy, which can 
then be released during star's evolution. Formation and 
structure of non-equilibrium neutron star crust was 
consi dered by many authors JVartanyan fc Ovakimova I 
11976'; 'Bisnovatvi-Kogan fc Chechetkin I Il979l; ISato I \l9m 



iHacnscl & Zdunik 1990, 2003; Gupt a et al. I l2007f) . Such 
, a state of the crust can be produced as a result of accretion 
' of the matter onto a neutron star in a close low-mass X-ray 
binary, where the original crust built of a catalyzed mat- 
. ter could actually be replaced by a crust with composition 
' strongly deviating from the nuclear equilibrium one. 

Heating due to non-equilibrium nuclear processes tak- 
ing place in the outer and inner crust of an accreting 
neutron star (deep crus tal heating) was calculated by 
IHaensel fc ZdunikT ([1990) who used a simple model of one- 
component plasma, and assumed that the outer layers of 
the matter produced in the X-ray burs ts consisted of pure 
^^Fe. Another simplification used by IHaensel fc Zduniki 
(|l990t ) consisted in assuming the ground-state - ground- 
state nuclear transitions due to the electron captures. 
Consequently, they maximized neutrino losses. Their cal- 
culated total deep crustal heating, produced mainly in the 
inner crust, was Qtn t ^1. 4 MeV per one accreted nucleon. 
[Hacnscl & Zdunik] ( 20031 ) recalculated deep crustal heat- 



ing for different initial composition of the outer layers, and 
obtained similar values of Qtot = 1-2 — 1.4 MeV/nucleon. 
Recently, heating of the outer crust of an acc reting neu- 
tron star was studied by 'Gupta et al. ' (2007) who went 
beyond a simple m o del o f Hacnscl fc Zdunik (19 90 ) and 
IHaensel fc Zdunik I (|2003D . Namely, iGupta et al. I (|2007t ) 
considered a multicomponent plasma, a reaction network 
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Up to now, the crustal heating used in modeling 
SXTs was described using the model of Haensel & Zdunik 
(1990) (hereafter referred to as HZ90), updated and gener- 
alized by Haensel & Zdunik (2003) (hereafter referred to as 
HZ03). In these models, the heat was produced during the 
active (accretion) episodes, when the outer layer of neu- 
tron star was sinking in the neutron star interior under the 
weight of accreted matter. The very outer layer was as- 
sumed to be composed of the ashes of the X-ray bursts in 
the active epoch. For simplicity, those ashes were assumed 
to be a one-component plasma {^^Fe in HZ90 and ^^Fe and 
i06pj HZ03). Under an increasing pressure, the composi- 
tion of a sinking matter element was changing in a sequence 
of nuclear reactions: electron capture, neutron emission and 
absorption, and finally, at densities exceeding 10^^ g cm"'^, 
also pycnonuclear fusion. As the nuclear processes were pro- 
ceeding off-equilibrium, they were accompanied by the heat 
deposition in the crustal matter. The crustal heating was 
mostly supplied by the pycnonuclear fusion processes in 
the inner crust a,t p — 10^^ — lO^'^ g cm^'^. This seem- 
ingly crucial role of pycnonuclear reactions is embarassing, 
because t heir rates are plagued b y huge uncertainties. As 
shown by lYakovlev et al71 (j2006[ ). the uncertainty in the 
calculated rate of pycnonuclear fusion of two ^*Ne nuclei, 
first pycnonuclear fusion in the inner crust as predicted 
by the HZ90 model, can be as large as seven (!) orders of 
magnitude. Therefore, there is a basic uncertainty about 
which pycnonuclear fusions do occur and at what densi- 
ties. Fortunately, as we show in the present paper, this 
uncertainty does not affect significantly the value of Qtot- 
However, it implies an ignorance concerning the distribu- 
tion of heat sources in the inner crust. Possible observa- 
tional constraints on the distribution of heat sources in the 
inner crust are discussed in Sect. [6] 

In the present paper we redo the calculations of the 
crustal heating taking into account uncertainties in the py- 
cnonuclear reaction rates as well as the suppression of the 
neutrino losses. We show that the uncertainties in the pyc- 
nonuclear reaction rates do not significantly affect the total 
heat release in the crust, while influencing of course the 
radial distribution of heat sources. We also show that the 
effect of excited states of daughter nuclei, while very im- 
portant in the outer crust, does not lead to a significant in- 
crease of total crustal heating. Calculations are performed 
assuming two different initial nuclides produced in the X- 
ray bursts. 

The plan of the paper is as follows. In Sect. [2] we briefly 
remind the scenario of the deep crustal heating, nuclear 
model used in simulations, and describe the relevant phys- 
ical processes acting in the crust of an accreting neutron 
star. Results of selected simulations of the nuclear evolu- 
tion of a matter element compressed from 10* gcm""^ to 
10^3.6 g^pg presented in the form of figures in Sect. 

[3] and tables in the Appendix. Total crustal heating is cal- 
culated in Sect.m In Sect. [S] we give an explanation of the 
weak model dependence of the total crustal heating (per one 
accreted nucleon). Our conclusions are presented in Sect.[Sl 
where we also suggest an observational testing of the actual 
radial distribution of heat sources, which could be helpful 
in putting constraints on the deep crustal heating models. 
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Fig. 1. (Color online) Integrated crustal heating per one 
accreted nucleon versus electron Fermi energy fj,e ■ The step- 
like curves were obtained using the HZ* model of the 
present paper, for two choices of an initial nuclide. The 
continuous curve was obtained bv lGupta et al. I (|2007[ ). 



2. Non-equilibrium nuclear processes 

2.1. Outer crust 

In what follows we briefly describe the nuclear evolution 
scenario of HZ90 and HZ03, w ith correction implied by the 
results of lGupta et al. I ()2007f ). Under conditions prevailing 
in accreting neutron-star crust at p > 10* g cm"'^ matter 
is strongly degenerate, and is relatively cold (T < 10* K), 
so that thermonuclear processes involving charged parti- 
cles can be assumed to be blocked by the Coulomb bar- 
rier. Consequently, the densities lower than the threshold 
for the pycnonuclear fusion (which is very uncertain, see 
lYakovlev et"ani2006l ppyc ~ lO^^-lO" g cm'^), the num- 
ber of nuclei in an element of matter does not change during 
the compression resulting from the increasing weight of ac- 
creted matter. Let us remind, that we assume that only one 
nuclear species {A, Z) is present at each pressure (one com- 
ponent plasma). Due to the nucleon pairing, stable nuclei 
in dense matter have even N = A — Z and Z (even-even 
nuclides). In the outer crust, in which free neutrons are 
absent, the electron captures which proceed in two steps. 



iA,Z) 
{A,Z-1) 



(AZ-l) 
{A,Z-2)- 



(1) 



lead to a systematic decrease of Z (and increase of = 
A — Z) with increasing density. The first capture in Eq. ([1]) 
proceeds in a quasi-equilibrium manner, with a negligible 
energy release. It produces an odd-odd nucleus, which is 
strongly unstable in dense medium and captures a second 
electron in an non-equilibrium manner, with energy release 
Qj, where j is the label of the non-equilibrium process. 
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2.2. Outer crust: c ompa rison of HZ* model with 
\GuDta et al. I (200^) 

In the original model of HZ90 electron captures were as- 
sumed to proceed from the ground state of the initial nu- 
cleus to the ground state of the daughter nucleus (GS-GS), 
and therefore most of the energy releas e was taken away by 
neutr inos (from 3/4 to 5/6 of Qj, see iHaensel &: Zdunik I 
l2003l for the discussion of this point) . Very recently, an ex- 
tensive, multicomponent self-consistent calculation of the 
nuclear evolution of an accreted matter eleme nt in the crust 
of a b ursting neutron star was carried out bv lGupta et al. I 
(|2007t ). As they have shown, electron captures to excited 
states (GS-ES) and subsequent de-excitation strongly re- 
duce neutrino losses, so that nearly all released energy heats 
the crust matter. We therefore modified the original HZ 
model by neglecting neutrino losses accompanying electron 
captures: this model will be denoted as HZ*. Consider a 
version HZ* model, corresponding to specific choice of the 
the initial mass number Ai, and denoted by a label (a). The 
integrated heat deposited in the crust in the layer with bot- 
tom density p can then be obtained as 



i(Pj<p) 



if 



(2) 



In Fig. [T] we compa re Q(ne), obta i ned fo r two HZ* mod- 
els, with results of iGupta et al. I (|2007( ). We replaced p 
variable by the electron Fermi energy / Ze to facilitate the 
comparison. The curve of lGupta et al. I (|2007D hes between 
two extreme HZ* models; it has a very large number of 
small jumps. These two features reflect the multicompo- 
nent structure o f crustal matt e r and large network of nu- 
clear reactions in lGupta et al. I (|2007f) . with 50 < < 110. 
Generally, viewing the simplicity of the HZ* model, its 
agreem ent with mult i comp onent, self consistent calcula- 
tions of iGupta et al. I ()2OO70 . is very good. 

2.3. Inner crust 

Above the neutron-drip point (p > pnd), electron captures 
trigger neutron emissions, 

{A,Z) + e- {A,Z-l) + iy,, 

{A,Z-l) + e- — > {A-k,Z -2) + kn + iy^ + Qj . (3) 

Due to the electron captures, the value of Z decreases with 
increasing density. In consequence, the Coulomb barrier 
prohibiting the nucleus- nucleus reaction lowers. This effect, 
combined with the decrease of the mean distance between 
the neighboring nuclei, and a simultaneous increase of en- 
ergy of the quantum zero-point vibrations around the nu- 
clear lattice sites, opens a possibility of the pycnonuclear 
reactions. The pycnonuclear fusion timescale Tpyc is a very 
sensitive function of Z. The chain of the reactions ^ leads 
to an abrupt decrease of Tpyc typically by 7 to 10 orders of 
magnitude. Pycnonuclear fusion switches-on as soon as Tpyc 
is smaller than the time of the travel of a piece of matter 
(due to the accretion) through the considered shell of mass 
Mshc\\{N, Z), Tacc = Afshou/-A/. The masses of the shells are 
of the order of 10~^ M©. As a result, in the inner crust the 
chain of reactions ([3]) in several cases is followed by the py- 
cnonuclear reaction on a timescale much shorter than Tacc- 
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Fig. 2. (Color online) Z and N of nuclei, versus mat- 
ter density in an accreting neutron-star crust, for different 
models of dense matter. Solid line: Ai = 56; dashed line: 
Ai = 106. Dash-dot line: continuation of Aj = 56 evolution 
but with pycnonuclear fusion blocked until Z = Zaun — 4. 
Each change of N and Z, which takes place at a constant 
pressure, is accompanied by a jump in density (see HZ90 
for detailed discussion of this point). Small steep segments 
connect the top and the bottom density of thin reaction 
shell. Arrows indicate positions of the neutron drip point. 

Denoting Z' — Z — 2, we have then 

{A, Z') + (A, Z') (2A, 2Z') + Q,- 1 , 

(2A, 2Z') — > {2A - k', 2Z') n + Q^- 2 , 



(4) 



where ". . ." correspond to some not specified chain of the 
electron captures accompanied by neutron emissions. The 
total heat deposition in matter, resulting from a chain of 
reactions involving a pycnonuclear fusion, is Qj = Qj^i + 
Qi.2 + Qi.3- In con trast to HZ03, and in accordance with 
iGupta et al~l (l2007h . we neglect the neutrino losses accom- 
panying the non equilibrium electron captures. 

Our model of atomic nuclei is described in detail in 
HZ90. Using our nuclear-evolution code, we evolved an el- 
ement of matter consisting initially of nuclei (Ai, Zi) im- 
mersed in an electron gas, from pi = 10^ g cm~'^ to 
p = Pf > 5 X 10^^ g cm~'^. Our results for the composi- 
tion and crustal heating are presented in the next section 
and in the Appendix. 

3. Results of simulations: reactions, heating, 
compositions 

The compositions of accreted neutron star crusts are shown 
in Fig. [2] and in tables in the Appendix. These results de- 
scribe crusts built of accreted and processed matter up to 
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Fig. 3. (Color online) Heat sources in the outer (upper 
panel) and inner (crust) for three HZ* models. Vertical 
lines, positioned at the density at the bottom of the re- 
action shell, represent the heat per one accreted nucleon. 
Labels as in Table 2. 
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Fig. 4. (Color online) Integrated heat released in the crust, 
Q{p) (per one accreted nucleon) versus p, assuming initial 
ashes of pure ^^Fe. Solid line: HZ* model of the present 
paper, with suppressed neutrino losses. Long dashes: GS- 
GS transitions in electron captures, with maximal neutrino 
losses. Dash-dot line: No neutrino losses, with pycnonuclear 
fusion blocked until Z — — 4. 



the density 5 x 10^^ g cm ^. At a constant accretion rate 
M = M_9 X 10"^ Mo/yr this wiU take ~ 10^ yr/M-g. 
During that time, a shell of X-ray burst ashes will be com- 
pressed from ~ 10^ g cm~'^ to ^ 10^^ g cm~'^. 

Two different compositions of X-ray bursts ashes at 
^ 10^ g cm"-^ were assumed. In the first case, Ai — 56, 
Zi = 26, like in HZ90. In the second scenario we followed 
HZ03, with A j = 1 06, to imitate nuclear ashes obtained by 
ISchatz et ahl (pOOlO . The value of Zi = 46 stems then from 
the condition of beta equilibrium at p = 10® g cm~^. The 
density dependence of Z and N = A—Z within the accreted 
crust is shown in Fig. O After pycnonuclear fusion region 
has been reached, both curves converged (as explained in 
HZ03, this results from Ai and Zi in two scenarios differing 
by a factor of about two). Suppressing pycnonuclear fusion 
in the Ai = 56 case until Z — Zmin = 4 yields the lowest 
curve. 

In Fig. [3] we show the heat deposited in the matter, 
per one accreted nucleon, in the thin shells in which non- 
equilibrium nuclear processes are taking place. Actually, 
reactions proceed at a constant pressure, and there is a 
density jump within a thin "reaction shell". The vertical 
lines whose height gives the heat deposited in matter are 
drawn at the density at the bottom of the reaction shell. 

One notices a specific dependence of the number of heat 
sources and the heating power of a single source on as- 
sumed Ai. In the case of Ai = 56 the number of sources 
is smaller, and their heat-per-nucleon values Qj are signif- 
icantly larger, than for Ai = 106. 



4. Integrated crustal heating 

The quantity Q^'^\p), for three specific models of compres- 
sional evolution, is plotted in Fig. [H In the all three cases, 
we set Ai = 56 and Zi = 26. For the first model, we ne- 
glect neutrino losses; its cumulated heat Q^^^ is always the 
highest. The second model is used to visualize the impor- 
tant of excited states of the daughter nuclei in the electron 
captures. For this model we assumed that the nuclear tran- 
sitions associated with electron captures are of the GS-GS 
type, which maximizes the neutrino losses. While the ef- 
fect is dramatic for p ^ 10^^ gcm~'^, it implies only a 20% 
underestimate of Q above 10^^'^ gcm^'^. But the most in- 
teresting is perhaps the effect of literally switching off of 
the pycnonuclear reactions, assumed in the third scenario. 
This was done by assuming that the pycnonuclear fusion is 
blocked till the nuclear charge goes down to Zmin = 4. And 
yet, for pb > 10^'^ gcm~^, Q'^^^ is very similar to that ob- 
tained in the first scenario, which was most advantageous 
as far as the crust heating was concerned. Namely, a miss- 
ing pycnonuclear heating at p ~ 10^^ g cm~^ is efficiently 
compensated by the electron captures accompanied by neu- 
tron emission within the density decade 10^^ — 10^^ g cm~^. 
The values of Q saturate above lO^'^'^ gcm~'^, where 80% 
of nucleons are in neutron gas phase. All in all, for three 
scenarios with Ai = 56, the total deep crustal heat release is 
(1.6 — 1.9) MeV/nucleon. For Ai = 106, numbers are shifted 
downward by about 0.4 MeV/nucleon. The summary of our 
results for the total heat release is given in Table [TJ 
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Table 1. Total crustal heating Qtot for Ai — 56 and Ai = 106. First line: HZ* model of the present paper. Second 
line: old HZ03 model with maximal neutrino losses. Third and fourth lines: results obtained when neutrino losses are 
suppressed and pycnonuclear fusion is blocked down to Zmin = 6 and Zmin = 4, respectively. 
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Fig. 5. (Color online) Baryon chemical potential fih{P) 
for different versions of blocking of pycnonuclear fusion. 
"Standard pycno" corresponds to HZ* model of the present 
paper, with = 56 . Three other curves correspond 
to pycnonuclear fusion suppressed until to Z = — 
8,6,4, respectively. The upper continous curves, which 

nearly coincide, are defined as Jl\!^\P) = ^J''b'\P) + 

Si(Pj <p) Q^"^- Dependence oiJl\^\p) on (a) is negligible, 

so that /Ib"^(^) ^ V'biP)- lowest smooth solid curve - 
cold catalyzed matter. 



5. Constancy of the total heat release 

Remarkably weak dependence of the total heat release 
in the crust, Qtot, on the nuclear history of an element 
of matter undergoing compression from ~ 10® g cm~^ to 
~ 10^^'^ gcm"'^, deserves an explanation. Let us there- 
fore study the most relevant thermodynamic quantity, the 
Gibbs free energy per nucleon, which is the baryon chem- 
ical potential, /ib. In the T = approximation, HhiP) = 
[S{P) + P]/ni,{P). Minimizing fih{P), at a fixed P, with re- 
spect to the independent thermodynamical variables {A, Z, 
mean free neutron density rT„, mean baryon density nb, 
size of the Wigner-Seitz cell, etc.), under the constraint of 
electric charge neutrality, rip = ne, we get the the ground 
state of the crust a given P. This state is called " cold cat- 



alyzed" matter, and its baryon chemical fJ-'^\P) is repre- 
sented in Fig. [5] by a solid black line. All other /x^,"'' (P) 
curves, showed in Fig. [51 display discountinuous drops due 
to non-equilibrium reactions included in a given evolution- 
ary model (a), and lie above the solid black one. This vi- 
sualizes the fact that non-catalyzed matter is a reservoir 
of energy, released in non-equilibrium processes that push 
down matter closer to the absolute ground state. In spite of 
dramatic differences between different fi^^ (P) in the region 
where the bulk of non-equilibrium reactions and heating oc- 
cur, P = (lO'^" — lO'^^'^) dyncm"^, these functions tend to 
the ground state one for P ^ 10'^^ dyncm"^. The general 

structure of different IJ,\^\P) is similar. At same P, their 
continuous segments have nearly the same slope. What dif- 
fers iii"\P) are discontinuous drops, by Qj"\ at reaction 



thresholds Pj"^ The functions fi'i^'iP) can therefore be 
expressed as 



(")/ 



A.(")(p)«Mb(n- E 

j(p,<p) 



(5) 



where /ib(^) a smooth function of P, independent of 
(a) (for explanation, see Fig. [5] and its caption). For P > 
2Q33.5 dyncm"^ the values of Qj are negligibly small, and 

all ^1^"^ (P) become quite close to the ground state line. This 



implies, that the sum 
precision, independent of (a) 



has to be, to a rather good 



6. Discussion and conclusion 

In the present paper we reconsidered a model 

of deep crustal heati ng, formulated originally in 

Haensel fc Zdunik I Il99d 120031) . Following iGupta et al.1 
20071 ) ■ we suppressed neutrino losses associated with 



electron captures. In this way we got the HZ* model of 
crustal heating that despite its simplicity (one compo- 
nent plasma) agrees nicely with re sults of self-c on sisten t 
multi-component calculations of iGupta et al. I (j2007[ ). 
Using the HZ* model we obtained the total crustal heating 
Qtot = 1-5 MeV and 1.9 MeV/nucleon, for the initial ashes 
consisting of ^°^Pd and ^^Fe, respectively. We studied the 
dependence of crustal heating on the location of the pyc- 
nonuclear fusion processes within the crust. It turned out 
that the total crustal heating (per one accreted nucleon) 
is quite insensitive to the depth at which pycnonuclear 
fusion occurs, with Qtot varying at most by 0.2 MeV, and 
we presented an explanation of this feature. 

Maximal neutrino losses, implied by the assumption 
of the ground-state - ground-state electron captures in 
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iHaensel fc Zdunik I (|1990f) and iHaensel fc Zdunik I (|2003f) . 

led to a severe u nderestimate of the heating of the outer 
crust (jGupta et a l. 2007). However, the outer crust gives 
a rather small contribution to the total crustal heating, 
Qtot, and the underestimate of Qtot is less than 25%. 
Composition of the initial ashes of the X-ray bursts can 
be more important for the v alue of Qtot, as shown already 
in (IHaensel fc Zdunik Il2003[) . 

The insensitivity of Qtot to the very uncertain rates of 
pycnoclear fusion means that this tremendous uncertainty, 
unlikely to be removed in spite of future theoretical efforts, 
does not affect the thermally equilibrated quiescent state 
of the SXTs. This may seem as good news. On the other 
hand, this means that the studies of SXTs in quiescence 
will not improve our knowledge of pycnonuclear fusion in 
dense plasma: this is bad news. 

Fortunately, the situation changes with the access to ob- 
servations of the thermal relaxation in SXTs after the accre- 
tion episode. This phenomenon cannot be observed in stan- 
dard SXTs that remain in accreting state for days - weeks. 
Thermal relaxation can be observed only in so called per- 
sistent SXTs, characterized by accretion states lasting for 
years - decades. Such thermal relaxation, called initial cool- 
ing, has been observe d in KS 1731-260 and MXB 1659-29 
(jCackett et al. II2006D . Let us focus on KS 1731-260. After 
12.5 years of accretion and associated crustal heating, the 
crust of KS 1731-260 has become significantly hotter than 
the neutron star core. After accretion stopped, the heat cu- 
mulated in the crust diffused over the star, and the stellar 
surface cooled. The cooling curve of KS 1731-260 toward th e 
quiescent state has been obtained bv lCackett et al. I (|2006D . 
This curve depends on the crust thermal conductivity, crust 
thickness, distribution of crustal heating sources, and on 
the neutrino cooling of ne utron star core (jRutledge et al. I 
[200l iCackett et al. l[2006h . 

There are two complementary aspects of observations 
of cooling of persistent SXTs. On the one hand, we need a 
theory to understand this phenomenon. On the other hand, 
observational data yield constraints on theoretical models. 
Very recently, simulations of thermal relaxa tion of KS 1731- 
260 w ere performed along these lines by IShternin et al. I 
(|2007l ). Increasing in number and precisions, observations 
of cooling curves in persistent SXTs will hopefully be a 
promising testing ground for the theories of deep crustal 
heating, including pycnonuclear fusion, and other physical 
processes in neutron stars. 
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Table A.l. Non-equilibrium processes in the crust of an accreting neutron stars assuming that the X-ray ashes consist 
of pure ^"^Pd. Pj and pj are pressure and density at which the reaction takes place. Ap/ pj is relative density jump 
connected with reaction, Qj is the heat deposited in the matter. Xn is the fraction of free neutrons among nucleons, and 
Pe is the electron chemical potential, both in the layer just above the reaction surface. 



p 


HJ 













( c\w\ cm ^ 








% 


(MeV) 


Q 9S'i X 1 n-^^ 

U'iJiJ'J i-\J 


3 51 7 X 1 0"'* 







4.4 


2.29 


22.9 




5 621 X 10°^ 


lOBp. 106 lyr n - , 





4.6 


6.34 


22.7 


9 ^79 X 1 Vp^ 


9 41 3 X 1 0'^'' 

^.t:J.O i-\J 


lOBlVf^ -^106 7^ _ Op- 1 2,, 





4.9 


10.43 


22.4 


8 581 X 1 n^* 


6 639 X 10^° 


lOSv^ 106 o _ n - 1 n 





5.1 


14.56 


22.2 


2.283 X 10^^ 


1.455 X 10^^ 


106g^ ^106 j^j. _ _^ IV, 





5.4 


18.73 


22.1 


5.025 X 10^^ 


2.774 X 10^^ 


lo^Kr ^'"^ Se - 26^ + 2v, 





5.7 


22.91 


22.0 


9.713 X 10^" 


4.811 X 10^^ 


loegg ^106 _ + 2v, 





6.1 


27.08 


22.0 


1 703 X 10'^" 


7 785 X 10'^ 


^""Ge Ni -1- 14n - 4e" + Av^ 


0.13 


13.2 


31.22 


110.8 


1.748 X 10^" 


8 989 X 10^^ 


^^Ni -^^^ Fe -1- 6n - 2e~ + 2ve 


0.19 


6.9 


31.33 


53.2 


1.924 X 10^" 


1.032 X 10'^ 


86pg ^80 Cj. _^ g„ _ + 2l/e 


0.25 


7.3 


31.02 


57.5 


2.135 X 10^" 


1.197 X 10'^ 


SOQj. ^74 _|_ g„ _ + 2z^e 


0.30 


7.7 


32.76 


62.1 


2.394 X lO^'^ 


1.403 X 10'^ 


74rp. ^68 Qg^ _^ _ 2e- + 2!/e 


0.36 


8.1 


33.57 


67.2 


2.720 X 10^° 


1.668 X 10'^ 


68ca ^''^ + 6n - 2e~ -I- 


0.42 


8.5 


34.45 


72.9 


3.145 X 10^° 


2.016 X 10'^ 


'^^Ar -^^"^ S -1- 6ri - 2e^ -1- 2j/e 


0.47 


9.0 


35.44 


79.2 


3.723 X 10^" 


2.488 X 10'^ 


56g ^50 Si _|_ gjj _ 2e- + 2v^ 


0.53 


9.4 


36.57 


86.0 


4.549 X 10^" 


3.153 X 10^^ 


^"Si Mg -t- 8n - 2e- ^ 2v,^ 


0.61 


8.8 


37.86 


94.5 


4.624 X 10^" 


3.472 X 10'^ 


*^Mg Ne -1- 6n - 2e" + 2i^e 














^''Ne Ne ^^^^ Ca 


0.66 


10.6 


37.13 


268.2 


5.584 X 10"'" 


4.399 X 10'^ 


'^Ca Ar + 6n - 2e" -I- 2i/e 


0.69 


4.8 


37.56 


31.6 


6.883 X 10^" 


5.355 X 10^^ 


66 ^60 s + 6ri - 2e" + 2i^e 


0.72 


4.7 


39.00 


33.5 


8.749 X 10^" 


6.655 X 10^^ 


60g ^54 Si _|_ g^ _ 2e~ + 2!/e 


0.75 


4.6 


40.33 


35.2 


1.157 X lO-"' 


8.487 X 10'^ 


^'^Si Mg + 8n - 2e- + 2i^e 














46]VIg +46 Mg 


0.79 


4.0 


41.84 


145.3 


1.234 X 10^' 


9.242 X 10'^ 


y^Cr Ti + 6n - 26" + 2v^ 


0.80 


2.0 


40.88 


11.4 


1.528 X 10^^ 


1.096 X 10'^ 


86rp. ^80 _^ g^ „ 2e" + 


0.82 


1.9 


42.04 


11.4 


1.933 X 10^^ 


1.317 X 10'^ 


^°Q,& -^''^ Ar + 6n - 2e" 2z/e 


0.83 


1.8 


43.31 


11.2 


2.510 X 10^^ 


1.609 X 10^^ 


74 Ar ^68 s 4. g„ _ 2e- + 2v^ 


0.85 


1.7 


44.71 


10.6 


3.363 X 10^" 


2.003 X 10'^ 


6»s si + 6n - 2e- +2i/e 














62si+62 Si^l^-* Ni 


0.86 


1.7 


46.26 


70.5 


4.588 X 10''' 


2.520 X 10'^ 


i^^Ni ^i^u Fe + 4n - 26" + 2i^e 


0.87 


0.8 


47.78 


3.0 


5.994 X 10^' 


3.044 X 10^3 


laope Cr -h 2n - 2e- + 2v,^ 


0.88 


0.9 


49.65 


2.7 


8.408 X 10^' 


3.844 X 10" 


ii8cr ->"6 Ti 2n - 26" + 2v^ 


0.88 


0.8 


52.28 


2.5 



Table A. 2. Nuclear processes and released heat in the inner crust, assuming initial ashes of pure ^°^Pd (i.e. as in Table 
lA.ip but suppressing pycnonuclear fusion until Z — Zmin — 4. Only the lines different than those in the Table lA.ll are 
presented. The network of reactions below pressure P = 4.624 x lO'^" dyn cm~^ is the same as in the Table ETT] 



(dyn cm ^) 


Pj 

(g cm-3) 


reactions 




% 


Pc 

(MeV) 


Qj 
(keV) 


4.624 X 10^" 


3.472 X 10'^ 


"^Mg Ne + 6n - 26" + 2i^e 


0.66 


9.9 


37.13 


76.0 


6.253 X 10^° 


4.745 X 10'^ 


36Ne „^30 + 6n - 2e- + 2v^ 


0.72 


9.9 


38.79 


84.2 


9.323 X 10^° 


6.937 X 10^2 


30o c + &n~2e- + 2i/, 


0.77 


9.4 


40.89 


91.8 


1.615 X 10^^ 


1.119 X 10'^ 


24q ^18 Bg _^ gj^ _ 2e- + 2Ue 


0.84 


9.2 


43.59 


94.4 


3.500 X 10-*' 


2.071 X 10''' 


"Be Li -h 3n - le~ + Ive 














i^Li -h^^ Li ^28 c + 2n 


0.88 


2.3 


46.78 


152.8 


6.339 X 10^' 


3.166 X 10" 


^"C^^' + In - le- + Ive 


0.89 


1.9 


48.98 


8.7 
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Table A. 3. Non-equilibrium processes in the crust of an accreting neutron stars assuming that the X-ray ashes consist of 
pure ^^Fe. Pj and pj are pressure and density at which the reaction takes place. Ap/pj is relative density jump connected 
with reaction, Qj is the heat deposited in the matter. X„ is the fraction of free neutrons among nucleons, and p,e is the 
electron chemical potential, both in the layer just above the reaction surface. 



p 


P 


Process 




Ap/p 


He 


q 


(dyn cm~^) 


(g cm-3) 










(keV) 


7.23 X 10^® 


1.49 X 10*^ 


56Fe ^^'^ Cr - 2e- + 2ue 





0.08 


4.08 


40.7 


9.57 X 10^^ 


1.11 X 10^° 


5^Cr Ti - 2e- + 2ve 





0.09 


8.18 


35.8 


1.15 X 10^^ 


7.85 X 10" 


56rpj __^56 _ + 2u, 





0.10 


15.64 


47.3 


4.75 X 10^*^ 


2.50 X 10" 


^^Ca. Ar - 2e' + 2v^ 





0.11 


22.48 


46.1 


1.36 X 10^° 


6.11 X 10" 


^^Av S + 4n - 2e~ -1- 2ve 





0.12 


29.38 


59.8 


1.980 X 10^° 


9.075 X 10" 


52g ^46 Si _|_ gjj _ 2e~ + 2Ue. 


0.07 


0.13 


32.27 


128.0 


9 V 1 n^o 


1 1 Qi V 1 n^^ 

_L . _L O X A ± U 


^''Si Mp- 4- fill — 9p^ 4^ 9?/ 


n 1 s 

U. -LO 


0.14 


32.22 




2 637 X 1 n'"*" 


1 455 X 1 0^^ 


™Me Np + 6t7 — 2p~ + 2u 














i>j e -p J-^ G — v^a 


u.oy 


u. r / 




ou / .y 


771 -v 1 n"^*^ 

Z. M i X iU 


1. ( DD X iU 


OcL — /\r -|- on - - ZG -j- ZUe 


U.40 


U.o 


Q/1 /1 7 


DO.O 


S 216 X 10^" 


2 1 "^4 X 10^^ 


62 A ^ ^56 Q 1 _ 4- 2z/ 


0.45 


0.09 


35.47 


71.6 


3.825 X 10^" 


2.634 X 10^^ 


56^^50 Si + 6n - 2e- + 2!/e 


0.50 


0.09 


36.59 


77.9 


4.699 X 10^" 


3.338 X 10^^ 


50si Mg 6n - 2e~ 2;/^ 


0.55 


0.09 


37.89 


84.6 


D.U-^iO X lU 


'^.o ( y X lu 


44TV/r ^36 at , o 9 ~ , r, 














SSlV-p ivr„ 72 














^^Ca -^^^ Ar + 6n - 2e' + 2i/e 


0.61 


0.14 


39.41 


308.8 


/ .ZOO X iU 


o.ooy X iu 


66 A „ ,60 c 1 0^— j_ 0,, 


U. ( u 


n n/i 


oy .ur 


9Q K 

zy .0 


Q 9"^^ X 1 n^" 

c/.ijtjo -i-\J 


7 041 X 1 0^^ 


6O0 ^54 C: 1 fi„ _ Op- -L 2?y 


0.73 


0.04 


40.34 


31.0 


l.ZZo X lU 


o.you X iu 


54q; ,48 T\ T 1 O-D— 1 0,, 

ol — ^ ~r D/l — Ze -|- ZUe 














4«Mg Mg Cr 














96 Cr ^^-^ Cr-h2n 


0.80 


0.04 


41.86 


135.1 


1.463 X 10^^ 


1.057 X 10^^ 


940- _^88 Ti + 6n - 2e- + 2iy^ 


0.81 


0.02 


41.99 


11.5 


1.816 X 10^^ 


1.254 X 10^^ 


ssrpj _^S2 Ca -f 6n - 2e" + 2ue 


0.82 


0.02 


43.18 


11.3 


2.304 X 10^^ 


1.506 X 10" 


®^Ca Ar 6n - 2e" + 2ue 


0.84 


0.02 


44.48 


10.9 


2.998 X 10^^ 


1.838 X 10" 


76 Ar -»™ S -h 6n - 2e- 4- 2i/e 


0.85 


0.02 


45.91 


10.0 


4.028 X 10^^ 


2.287 X 10" 


70s Si 6n - 2e- + 2v^ 














64si+64 gj^l28 














128Ni^l26 Ni-^2n 


0.87 


0.01 


47.48 


67.3 


5.278 X 10^^ 


2.784 X 10" 


i26Ni _^i24 pe + 2n - 2e- + 2ve 


0.88 


0.01 


48.50 


2.5 


7.311 X 10^^ 


3.493 X 10" 


i24pg ^122 Qj. _^ 2n - 2e- 2;/^ 


0.89 


0.01 


51.05 


2.4 



Table A. 4. Nuclear processes and released heat in the inner crust, assuming initial ashes of pure ^^Fe (i.e. as in Table 
EU but suppressing pycnonuclear fusion until Z = Zmin = 4. Only the lines different than those in the Table IA.3I are 
presented. The network of reactions below pressure P — 2.637 x lO'^" dyn cm^^ is the same as in the Table ETSl 



p 

(dyn cm~^) 


P 

(g cm-3) 


reactions 








Ap/p 

% 


Pc 

(MeV) 


1 

(keV) 


2.637 X 10™ 


1.455 X 10'^ 


*"Mg Ne + 6n 


-2e 


" + 2Ue 


0.40 


16.0 


34.45 


161.4 


3.227 X 10™ 


1.961 X 10^^ 


^'^Ne ^2** 6n - 


- 2e' 


+ 2Ve 


0.50 


17.9 


35.71 


183.4 


4.254 X 10™ 


2.831 X 10^^ 


28q ^22 c + 6n- 


2e" -1- 2!^e 


0.61 


19.4 


37.49 


209.8 


6.392 X 10™ 


4.541 X 10" 


^^C-^" Bc + 6ri- 


2e- 


+ 2ue 


0.72 


20.0 


39.89 


238.2 


1.216 X 10^' 


8.617 X 10'^ 


"Be->" Li + 4n- 

13Li+13 Li ^26 c 


- le" 


+ lUe 














26C ^24 c + 2n 






0.79 


7.5 


43.18 


331.2 


1.727 X 10^^ 


1.186 X 10" 


24c ^" Be + 6n - 


2e- 


+ 2ue 


0.85 


7.4 


43.64 


84.1 


3.798 X lO''' 


2.200 X 10'^ 


^«Be^^' Li-f In- 
"U+^^ Li ^^^C 


- le- 


+ IVe 














34^ ^30 C _^ 4^ 






0.88 


0.028 


46.86 


142.0 





